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Abstract 

The transitions in disordered substances are discussed briefly: liquid-liquid 
phase transitions, liquid-glass transition and the transformations of one amorphous 
form to another amorphous form of the same substances. A description of these 
transitions in terms of many-particle conditional distribution functions is proposed. 
The concept of a hidden long range order is proposed, which is connected with 
the broken symmetry of higher order distribution functions. The appearance of 
frustration in simple supercooled Lennard-Jones liquid is demonstrated. 



It is well known for a long time that there exist sharp phase transitions between differ- 
ent aggregate states and sharp polymorphic phase transitions between different crystalline 
phases of the same substances. During last two decades a lot of experimental data was 
obtained on complicated phase diagrams of liquids and amorphous solids, too. Some of 
these results were presented on the first international conference that took place here in 
Russia in 2001 PQ and that can be considered as a formal claim of a new direction of 
physical investigations - transformations in disordered substances: liquid-liquid transi- 
tions and transformations of one amorphous form to another amorphous form of the same 
substances. The liquid-glass transition, although has longer history, has to be considered 
in the same context. It is now firmly established by different experimental techniques 
that sharp liquid-liquid transitions under pressure, formally similar to first-order phase 
transitions, exist as well as reversible transformations between amorphous states involving 
changes in local order structures and density. Usually a crystal melts with a conservation 
of the short-range order (SRO) structure type or into denser liquid with SRO structure 
similar to that of high pressure crystalline phase. 

It should be emphasized that the transitions in liquids are true phase transitions mainly 
determined by thermodynamic relationships, whereas the transitions in amorphous solids 
take place far away from equilibrium and are governed by the corresponding kinetics. 

A useful microscopic theory of these transitions is not developed yet and only empir- 
ical models and computer simulations have been used in practice to date. For example, 
interesting results were obtained by Stanley through molecular dynamics study basing 
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on the taking into account the hydrogen bonding in supercooled water [2|. It is by the 
demonstration of a simple analytic way of obtaining Stanley results that we begin the 
presentation of our own results on this subject (see, e.g., [TJ El IH El E] and references 
therein). 

From the intuitive point of view liquid-liquid phase transition between low density and 
high density phases may be related to the competition between expanded and compact 
structures. This suggests that the potential should have two equilibrium positions. The 
most obvious form of such potential is: 



$(r) 



This two-well potential may be considered as a model for the water potential 

To investigate the possibility of the existence of the second critical point in this case we 
developed the mean-field (van der Waals-like) theory. Using the well-known Bogoliubov 
inequality for the free energy we can write F < F HS + < U — Uhs >hs- Here F HS is the 
free energy of the system of hard spheres with diameter a, and we consider the attractive 
part as a perturbation. Here U = ~ ®( r ij) an d Uhs — \ 5_/«<y ®Hs( r ij)- The average 
over the hard sphere potential has the form 

POO 

<U-U HS >hs= 2npN / <£ atr (r)g H s(r)r 2 dr, (2) 

Jo 

where $ a4r (r) = $(r) — $#s(r), gns{ r ) is the radial distribution function of the hard 
sphere system which we take in the Percus-Yevick approximation [7] and for Fhs we use 
the approximate Carnahan-Starling equation [Hj 

FHS 31nA-l + lnp+ 'r 3 l (3) 
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Here A = h/(2iTmk B T) 1 / 2 . 

The equation of state is given by P = p 2 d(F/N)/dp. In Fig. 1 the two families 
of isotherms are shown for the temperatures close to two critical points (3\ = £i/&bT c i 
and fa = £i/&b7c2- Fig. 1 shows that at low temperatures (/3 = ei/ksT > flx) there 
are two van der Waals-like loops in the equation of state which correspond to two fluid- 
fluid transitions. In the temperature region /3 2 < (3 < (3% there is only one loop which 
corresponds to the well known gas-liquid transition, /3 2 being the gas-liquid critical point 
temperature and fix - the liquid-liquid critical temperature. 

In this example, as at the ordinary critical point, no symmetry of the correlation 
functions is broken at the transition. The order parameter is the difference of densities 
of high and low density phases Ap = pn — p i2 . However, it is interesting to describe the 
change of the local structure and the cluster symmetry at the transition. 

To describe different kinds of space symmetry breaking we use the formalism of clas- 
sical many particle conditional distribution functions 

p ( r i_o r o\ _ -fl»+i(ri,rg,...,rg) 
F ' +l(ri|r "" r * ) F,(r» r») ' 
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Here F s (r 1; r a ) is the usual s-particle distribution function 0. The functions F s+1 (r 1 |r ( j ) ...r°) 
satisfy the equation 



pF s+ i(ri|r t | ) ...r ^ 



exp|-/?^$(r 1 -r°) + ^^- / S k+1 (r u ...,r k+1 ) 
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k=l k>l 

F s+1 (r|r?...r°)...F s+1 (r fe+1 |r?...r°)dr 2 ...dr fe+1 }. (4) 



Here z is activity, p is the mean number density, S k +i(ri, r^+i) is the irreducible cluster 
sum of Mayer functions connecting (at least doubly) k + 1 particles, (3 = 1/HbT. 

The simplest case is the symmetry breaking of the one-particle function. In the solid 
phase the local density, proportional to the one-particle distribution function, has the 
symmetry of a crystal lattice and can be expanded in a Fourier series in reciprocal lattice 
vectors G: 

p(r) = ]Tp G e* Gr , (5) 

G 

where the Fourier coefficients pc are the order parameters for the transition. 

The Taylor expansion of the corresponding free energy functional around the liquid 
can be written in the following form: 

PAF= ! drg(r)\n^-J2y [ c (n) (rx, r k )Ag(r 1 )...Ag(r k )dr 1 ...dr k , (6) 

J ^° k>2 ' J 

where 

Ag(r) = g(r) - Qi 

is the local density difference between solid and liquid phases. This is the base of the 
density functional theory of freezing (DFT) j^j. In the frame of this approach tens of 
melting curves were calculated (see, e.g., the reviews [TT]). The full system of equations to 
be solved in DFT contains the nonlinear integral equation for the function p(r), obtained 
as the extremum condition for the free energy and the equilibrium conditions for the 
chemical potential and the pressure written in terms of the same functions as in To 
proceed constructively in the frame of DFT we have to choose an actual form of the free 
energy functional - a kind of closure or truncating - and we must make an ansatz for the 
average density of the crystal. The importance of such an ansatz follows from the fact 
that we are dealing with a theory which is equivalent to Gibbs distribution and one has 
to break symmetry following the Bogoliubov concept of quasiaverages [T2] . 

Now let us consider a state of matter which is characterized by the uniform local 
density, but the broken symmetry of the two-particle distribution function. Such type of 
order is called the bond orientations order (BOO), where "bond" is the vector joining a 
particle with its nearest neighbor. This kind of order is well known in theories of two- 
dimensional melting (hexatic phase) Near the transition to anisotropic liquid state 
we have: 

F 2 (r 1 |r?)=^(|r 1 -r?|)(l + /(r 1 -r?)), (7) 

where f(r\ — r°) has the symmetry of the local environment of the particle at r° and 
may be written in the form f{r\ — r®) = f(\ri — r°|,fi), Q determines the direction of 
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the vector r x — rj. In the case of three dimensions function f(r, Q) may be expanded in 
a series in spherical harmonics: 

oo I 

/(r,fl) = EE WO^)- (8) 

i=0 m=-Z 

The microscopic equations for the order parameters fi m { r ) can be obtained from the 
main equation (Jlj). The linearized equation determines the instability of isotropic liquid 
against the formation of the state with BOO and has the form 

Air f 

fim(r) - ^— j J r,(r, r')g(r')f lm (r'y 2 dr' = 0. (9) 
Here Ti(r,r') correspond to the isotropic liquid when 

r(r 1; r?,r 2 ) = £ -^-j | 5 fc+1 (n...r fc+1 ) x 

x ( 7 (|r 3 -r?|)--^(|r fe+1 -r°|)dr 3 ---rfr fe+1 . (10) 

reduces to 

r(r 1 ,r?,r 2 ) = r(r,r / ,0), (11) 



oo . I 

r(r,r',e)^J2^h r ^ r ') £ ^(^i)^(^), (12) 

1=0 l=-m 

The angles f2 x and Q2 determine the directions of the vectors r and r' and r = |ri — r°|, r' = 
|r2 — rJI, ^ is the angle between vectors r and r'. It should be notice that the correlation 
length of the orientational fluctuations £z >m — > 00 when approaching the instability line 
given by Eq. 0. 

To describe liquid-liquid and liquid-glass transitions we must consider isotropic case 
with rotationally invariant two-particle distribution function. A possible description of 
these cases can be given in terms of broken symmetry of higher order distribution func- 
tions. At high temperature the nearest neighbors of a molecule can take different relative 
positions and there is no SRO. At lower temperature SRO appears which can be of dif- 
ferent kinds at different densities. The rotation and the translation of the clusters of 
preferred symmetry give rise to the fact that one-particle and two-particle distribution 
functions remain isotropic. If a kind of BOO appears the clusters are oriented in similar 
way and the two-particle distribution function becomes to be anisotropic (as in 2D hexatic 
phase). However, we can imagine another situation - freezing of the symmetry axes of the 
clusters in different position. The isotropic phase can be considered as analogous to the 
paramagnetic phase (of cluster symmetry axes), the BOO phase - to the ferromagnetic 
phase, and the mentioned freezed phase - to a spin glass phase. 

Let us consider for simplicity a 2D system. In the vicinity of the transition one can 
write (in the superposition approximation for the liquid) 

F 8 (n|rS, rjj) = 0(|n - r?|M|n - r°|)(l + / 3 ( ri |r?, rjj) (13) 
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In 2D case / 3 (r!|r°,r° 




(14) 



where Ro = r° — r°, Ri = ri — r®, R 2 = r 2 — and 0o is the angle of the vector R with 
the z axis, 61 - the angle between Ri and Ro and © 2 - the angle between R 2 and Ro. 
The linearization of (@J) for s = 2 gives: 



f 3 (RoAo;Ri,<3>i)= / T'(R , ( j )0 ;r 2 ;R 1 ,e 1 )f 3 (Ro,cj ) o;R2,e2)g(\^2-'Ro\)g(R2)dr 2 , (15) 



There are two kinds of angles entering the equations and two kinds of order parameters, 
consequently. One angle (0 O ) h xes the position of one pair of particles of the cluster, 
and the other (0,) - the position of the third particle in the coordinate frame defined by 
0o- The order parameter connected with Bj is the generalization of intracluster hexatic 
parameter for the case of different coordinate frames. The order parameter connected 
with 0o is an analogue of magnetic moment and in glass-like phase one can consider an 
Edwards- Anderson parameter < cos0o(t) cos0o(O) >. In such a way we come to the 
concept of a "conditional" or "hidden" long range order: if we consider two pairs of 
particles at infinite distance from one another then there exists a preferable possibility 
for the relative position of the third particle near each pair. The directions of the bonds 
in the pairs of particles themselves are subjects to spin-glass-like order. In 3D case the 
rotation of clusters is given by rotation matrices D\™ (uj 0i ) so that we obtain a kind of 
orientational multipole glass for the clusters. If the intracluster ordering is established 
then we can consider the system of clusters. The orientational state of this system is 
defined by the intercluster interaction for different values of temperature an pressure. 

Now let us consider this later situation when the intracluster symmetry is fixed and let 
us try to estimate the intercluster orientational interaction. If the intercluster interaction 
had the same sign for all cluster sizes (or all clusters had the same size) one would get the 
state with simple BOO. However, because of the difference in cluster sizes the orientational 
interaction for some harmonics may change sign as a function of the cluster size (see Fig. 
2). In this case the low temperature state should be amorphous for some harmonics. So 
the difference of the orientational interaction of the clusters for different cluster sizes may 
be considered as the reason of some kind of frustration in simple liquids. It should be 
emphasized that the form of the corresponding component of the orientational interaction 
is the intrinsic statistical property of the liquid and does not depend on the timescale of 
the fluctuations in size and symmetry of the cluster. There is no real quenched disorder 
in the system but only an analog of it which may be treated in a formally same way as the 
quenched disorder in spin glasses. To analyze qualitatively the orientational freezing in 
the system we introduce simple lattice model which takes into account the interaction only 
between clusters with definite symmetry. The model gives the possibility to conclude what 
harmonics freeze first and what local symmetry prevails immediately below the transition. 
Let us now describe our results in more detail. 



where 




fc>i 



x g(\r 3 -r° 2 \)...g(\r k+1 -r° l \)g{\r k+1 -r 2 \)dr 3 ...dr k+1 . (16) 
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Our starting point is the expression for the free energy of the system as a functional 
of a pair distribution function g 2 (r,, r ) which has the form jlj: 



F/k B T = J drdr o /?02(r,r o ) [in (\ 3 pg 2 (r, r )) - l] - 
~E (n + l)\ / S n+i( r i-- r n+i)g2(r 1 ,r )---g 2 (r n+1 ,r ) x 

n 

xdrt ■ • ■ dr n+1 dr -J $(r - r )pg 2 {r, r )drdr , (17) 

where the term with logarithm corresponds to the entropy and the other terms — to the 
interaction energy. Here $(r — r ) - interparticle potential (for Lennard- Jones potential, 
$ (r) = Ae((a/r) 12 - (a/rf)), X = h/(2nmk B T) 1 / 2 . 

We can estimate the change of the energy due to (|7Jl. Omitting the entropy term in 
EqlTTT we have up to the second order in 5g(r, r ). 



AF/k b T = -i J T(r 1 ,r ,r 2 )Sg(v 1 ,v )Sg(r 2 ,r )dr 1 dr 2 . 



Using the approximation for the radial distribution function g(r) = p 1 (n s /4:7rr 2 )5(r — r s ) 
and the Eq.(jHJ) we obtain: 

AF(r s )/k B T = - l -p^ (^) 2 jt^ri(r s ,r s ) ^ f Y lm {^i) x 

1=0 m=-l 
1 X ' 

xY; m (Q 2 )f(r s ,Q 1 )f(r s ,Q 2 )dQ 1 dQ 2 = --^ J,(r s ) ^ |/ /m | 2 . (19) 

Z=0 m=-l 

Here Ji(r s ) = p~ 2 2^-(|^) 2 r/(r s , r s ), n s is the number of nearest neighbors of a particle 
and r s is the size of the cluster, which is of the order of the first coordination shell size. 

The function AF(r s ) may be interpreted as the mean-field orientational interaction 
energy of the system of clusters having the size r s . To get the full energy of the system 
one should integrate (fE?j) over the probability of finding the cluster with the size r s which 
is given by the function r 2 g(r s ) in the vicinity of the first maximum. 

Using the approximations of |1] for T(ri, r°, we obtain the estimation for Ji(r s ) as 
a function of r s . Fig.2 represents Ji(r s ) for I = 4 and 6 along with r 2 s g{r s ) in the vicinity 
of its first peak. It is seen that J/(r s ) changes sign. This result enables us to suppose 
that there is a kind of frustration (which is analogous to that in spin glasses) appearing 
as a result of variations in the sizes of clusters according to g(r) and that it is possible 
to study the transition in the system of interacting clusters on the base of simple model 
lattice Hamiltonian: 

H = - a £ E 4 E u ^muL(^)- (20) 

<i^=j> l=0 m=-l 

The functions L^ m (flj) are the lattice harmonics for the point groups corresponding to 
the cluster symmetry. This Hamiltonian describes correctly BOO of clusters. In this case 
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Figure 2: Ji(r s ) for / = 4 and 6 along with r 2 s g(r s ) in the vicinity of the first peak as 
functions of r s for dimensionless temperature ksT/e = 1.0 and density pa 3 = 1.2. 



the energy calculated from Eq. (f2*U|) in the mean-field approximation (taking into account 
that (Ui m (£li)) = fim) coincides with the intercluster energy (fTUj) under appropriate choice 
of J\y We will use the Hamiltonian (|20jl to study the system of interacting clusters with 
various sizes. 

To simplify the problem we neglect in Hamiltonian f!2()l) all the terms except ones cor- 
responding to the unit representation of the point group. Furthermore, we consider only 
the cases I = 4 and I = 6 which represent the cases of cubic and icosahedral symmetries. 
This Ising-like model may be called a "minimal" model: 

H = - 1 -J2^Au r (21) 

Functions U = U((p, 9) are the combinations of spherical harmonics. We will consider 
separately symmetries of "simple" cube (I — 4, m — 0, ±4), cube (I — 6, m = 0, ±4) and 
icosahedron (I = 6,m = 0, ±5) correspondingly [TH lib]. For example, for I = 4 one has: 

U = C%, 9) = \l — { Y m {^ 9) + X /^ (Y u {<p, 9) + r 44 (-y?, 9)) } (22) 

The interactions are chosen in such a way that the MF approximation gives exact 
solution (infinite-range interactions). It is easily seen that in the minimal model (|21|) 
without disorder in the framework of the MF approximation there is a first order phase 
transition to the state with BOO (compare to [T^ ITo]). From Fig. 2 it is clear that, as the 
first qualitative step, J|- may be chosen as random interactions with Gaussian probability 
distribution 

(Jjj — Jo) 
v^j [ 2 J 2 




P(Jij) = 7= — ^ exp 



(23) 



S 



where J = J/y/N , J = Jq/N can be related to the microscopic parameters. We 
approximate r 2 g(r) by a gaussian exponential near the position of the first maximum r . 
So r 2 g(r) ~ exp[— (r — r ) 2 /2a}. The approximation for the functions Ti is then linear: 
Ti ?z a + (3{r — r ). That is: Jq = a, J = (3^fa. 

The free energy of the system can be obtained using replica approach (see, e.g., [*"3~j). 
In the replica-symmetric (RS) approximation we have [T7j : 

kT ) 2 4 4 




(-^lnTr[exp (#)]}, (24) 



where the trace in this case is defined as follows: Tr(. . .) = f Q * dip dcos(6)(. . .). Here 
t = J/ksT and 
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The order parameters are: m is the regular order parameter (an analog of magnetic 
moment in spin glasses), q is the glass order parameter and p is an auxiliary order pa- 
rameter. The extremum conditions for the free energy (J24)) give the following equations 
for these order parameters: 



m={U), p=(U 2 ), q = (U) 2 , (25) 



where (. . .) = Tr(. . . e e )/Tre e and (. . .) = jT OQ '^ e ~ z //2 [---]- We find from these equations 
the temperature dependence of the order parameters. The RS solution is stable unless 
the replicon mode energy A rep i is nonzero [T"""l IT7] . For our model we have 



A repl = l-t 2 ((f/ 2 » 2 , (26) 

where ((. . .)) denotes the irreducible correlator. We find the temperature T A _ T that 
corresponds to A rep i = 0. To obtain the actual glass transition temperature one has to 
study the dynamics of the system. In this paper we limit ourselves by the static approach. 
As is usually believed 01201 and correctly shown in [21] the dynamical T g can be obtained 
in the frame of the static approach as the temperature T m of the marginal instability of 
the one-step RS breaking solution. We have calculated T m and found that within the 
accuracy of calculations T m and T A _ T coincide. We expect that as in spin glasses below 
T A _ T the liquid dynamics is characterized by long relaxation times and other phenomena 
characteristic to glass transitions. So there is the glass transition in the simple cube case 
with T A _ T ~ 0.39; in the other cases, icosahedron and cube, there is no glass transition but 
just a first order transition to BOO state at temperatures about 0.45, 0.42 correspondingly 
at per 3 = 0.973. The last two temperatures of BOO transitions are in agreement with the 
results of molecular dynamics simulations of Ref.[22j. It should be noted that all these 
temperatures are well below the melting temperature T = 0.703 at this density [2*2*1 H5| . 

The work was supported in part by the Russian Foundation for Basic Research (Grant 
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